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Abstract 

Transient  analytical  solutions  are  presented  for  the  overpotential  and  the  voltage  for  a  porous  electrode  that  includes  both  double-layer 
charging  and  a  faradaic  reaction.  In  addition,  a  simplified  dynamic  model  is  developed  for  the  same  process  based  on  a  second  order, 
three-parameter  polynomial  approximate  model.  The  effects  of  the  parameters  such  as  dimensionless  exchange  current  density,  conductivity 
ratio  and  applied  current  density  on  the  voltage  and  overpotential  distribution  are  presented.  Also,  using  the  exact  transient  solution  an 
expression  for  the  dimensionless  interfacial  current  density  is  derived  and  the  effects  of  the  parameters  mentioned  above  are  presented. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Capacitors  have  received  considerable  attention  recently 
because  of  their  high  power  densities  [1].  Originally,  capac¬ 
itors  served  mainly  as  backup  power  sources  for  memories, 
microcomputers  and  clocks.  Later  they  were  used  as  main 
power  sources  in  applications  that  required  short  pulses  of 
large  current  for  several  seconds  duration.  Currently,  they 
are  employed  in  hybrid  configurations  with  batteries  and 
fuel  cells.  Capacitors  can  be  broadly  classified  [2]  as  electro¬ 
chemical  double-layer  capacitors  or  pseudocapacitors  (i.e., 
an  electrochemical  capacitor  with  a  faradaic  reaction). 

The  performance  of  a  capacitor  depends  on  many  factors 
such  as  the  material  properties,  cell  design  and  operating 
conditions.  To  fabricate  capacitors  for  a  particular  applica¬ 
tion,  it  is  necessary  to  analyze  the  performance  and  optimize 
the  system  parameters.  Mathematical  modeling  has  proven 
to  be  very  useful  for  this  purpose  [3].  Both  circuit  analog 
models  and  macroscopic  models  [2-5]  have  been  used  to 
predict  the  behavior  of  capacitors.  Circuit  analog  models 
do  not  capture  the  physics.  While  macroscopic  models  cap¬ 
ture  the  physics,  they  are  often  time  consuming  to  solve. 
Hence  simplified  models  without  compromising  accuracy 
are  highly  desirable.  In  this  paper,  an  approximate  model 
for  a  pseudocapacitor  is  presented. 
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A  pseudocapacitor  cell  can  consist  of  two  porous  elec¬ 
trodes  with  a  separator  between  them  as  illustrated  in  Fig.  1. 
The  system  is  filled  with  electrolyte  throughout.  In  this  pa¬ 
per,  we  present  a  model  for  a  single  porous  electrode.  The 
objectives  of  this  work  are  three-fold:  (1)  Develop  an  analyt¬ 
ical  solution  for  the  transient  behavior  of  a  porous  electrode 
with  both  double-layer  charging  and  a  faradaic  reaction  with 
linear  kinetics.  (2)  Develop  a  simplified  dynamic  model 
to  replace  the  transient  solution  obtained  in  (1).  (3)  Study 
the  effect  of  the  various  parameters  such  as  the  exchange 
current  density,  conductivity  ratio  and  applied  current  den¬ 
sity  on  the  developed  model.  The  interfacial  current  density 
distribution  across  the  electrode  is  a  measure  of  the  active 
material  utilization.  Based  on  the  interfacial  current  den¬ 
sity  distribution,  the  performance  of  the  capacitors  can  be 
optimized.  Here,  an  expression  for  the  interfacial  current 
density  is  derived  using  the  transient  analytical  solution. 
Using  this  expression  the  interfacial  current  distribution  is 
determined  as  a  function  of  the  system  parameters. 


2.  Model  description 

The  following  assumptions  are  made:  (1)  Porous  elec¬ 
trode  theory  in  one  dimension  is  applicable.  [6].  (2)  No 
concentration  gradients  exist  inside  the  electrode.  (3)  Both 
double-layer  charging  and  a  linear  faradaic  reaction  occur. 
(4)  The  material  and  kinetic  properties  (< a ,  <r,  /c,  a a,  ac 
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Nomenclature 

a  specific  interfacial  area  (cm-1) 

a(j)  time-dependent  parameter  of  the 

simplified  dynamic  model  (see  Eq.  (25)) 
b( r)  time-dependent  parameter  of  the 

dynamic  model  (see  Eq.  (20)) 
c( r)  time-dependent  parameter  of  the  dynamic 

model  (see  Eq.  (21)) 

Cdi  double-layer  capacitance  (F/cm2) 

F  Faraday’s  constant,  96487  C/eq. 

/o  exchange  current  density  (A/cm2) 

i\  matrix  phase  current  density  (A/cm2) 

22  solution  phase  current  density  (A/cm2) 

I  total  cell  current  density  (A/cm2) 

jn  interfacial  current  density  (A/cm2) 

j*  dimensionless  interfacial  current  density 

(see  Eq.  (43)) 

jn, f  faradaic  interfacial  current  density  (A/cm2) 

L  thickness  of  the  porous  electrode  (cm) 

n  number  of  electrons  transferred  in  the  reaction 

R  universal  gas  constant,  8.313  J/molK 

t  time  (s) 

T  temperature  (K) 

V*  dimensionless  voltage  drop  across  the 

porous  electrode 
v  distance  (cm) 

X  dimensionless  distance,  x/L 

Greek  symbols 

aa,  ac  anodic  and  cathodic  transfer  coefficients 
respectively  (aa  +  ac  =  n) 
ft  ratio  of  the  solution  phase  and  matrix 

phase  conductivities,  k/o 

8  dimensionless  current  density,  I(FL/kRT) 

0  overpotential  (0 1  —  02)  (V) 

77*  dimensionless  overpotential,  r)F/RT 

77*vg  average  dimensionless  overpotential 

(see  Eq.  (24)) 

k  solution  phase  conductivity  (S/cm) 

v2  dimensionless  exchange  current  density 

[aio(aa  +  ac)FL2/RT](l/cr  +  1/k) 
o  matrix  phase  conductivity  (S/cm) 

r  dimensionless  time  t/(aC&(\ /k  +  1  /cr)L2) 

01  solid  phase  potential  (V) 

02  solution  phase  potential  (V) 

0*  dimensionless  solid  phase  potential,  0i  F/RT 


and  /o )  are  assumed  to  be  constants.  (5)  The  double-layer 
capacitance  (Cdi)  is  taken  to  be  a  constant.  (6)  The  open 
circuit  potential  is  set  equal  to  zero.  With  no  concentra¬ 
tion  gradients,  the  matrix  phase  current  density,  i\  and  so¬ 
lution  phase  current  density,  22  are  given  by  Ohm’s  law 
[6]: 


h 


901 

=  —<7 - 

dx 


h 


902 

=  —K - 

dx 


(1) 

(2) 


where  o  and  k  are  the  matrix  phase  and  solution  phase  con¬ 
ductivities,  respectively.  The  total  current  density  I  is  the 
sum  of  the  matrix  and  solution  phase  current  densities: 


21  +  22  =  / 


where  I  is  the  applied  current  density  under  galvanostatic 
conditions  (note  that  I  is  defined  to  be  positive  when  leaving 
the  current  collector  at  v  =  L,  see  Fig.  1).  The  current 
transferred  from  the  matrix  phase  to  the  solution  phase  is 
expressed  in  terms  of  the  interfacial  current  density,  jn  [6] : 


9/i 

dx 


where  a  is  the  surface  area  per  unit  volume  of  the  porous 
electrode.  The  interfacial  current  density  (jn)  is  the  sum  of 
the  double-layer  charging  current  density  and  the  faradaic 
current  density: 


jn  —  G  di 


9(01  ~  02) 

dt 


+  jn,  f 


where  Cdi  is  the  double-layer  capacitance  and  the  faradaic 
current  density  (jn  f )  is  given  by  the  linear  kinetic  expression: 


ipFjfx a  +  Q?c) 

RT 


where  aa  +  ac  =  n  [6]  and  the  open  circuit  potential  has 
been  set  equal  to  zero.  Substituting  Eq.  (6)  in  Eq.  (5)  yields 


9(01-02)  ,  io(oia  +  ac)F 

Jn  =  Cdi - - - 1 - — - (01  -  02) 


dt 


RT 


Let  the  overpotential  be  given  by  27  =  0i  —  02.  The  initial 
and  boundary  conditions  for  the  overpotential  are  given  by 

[5]: 


01  =  0  and  02  =  0  ^  7/  =  0  at  /  =  0  for  0  <  v  <  L 

(8) 


Positive  Negative 

Electrode  Separator  Electrode 


► 


Matrix  Phase 


*• 


Solution  Phase 


Current  Collector 


►  I 


x  =  0  x  =  L 


Fig.  1.  Schematic  of  an  electrochemical  capacitor. 
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i  i=0  and 


dr)  I 

i2  =  /  =>-  —  =  -  at  .y  =  0  for  t  >  0 
dx  K 


dr)  I 

i\  =  I  and  i2  =  0  =>  —  = - at  x 

dx  a 

=  L  and  for  t  >  0 


(10) 


2.2.  Simplified  dynamic  model 

An  approximate  solution  for  the  same  model  (Eq.  (11) 
subject  to  Eqs.  (14)— (16)  can  be  found  by  assuming  that  the 
dimensionless  overpotential  is  a  polynomial  in  X  [8]: 

r)*{X,  r)  =a(r)  +  b(t)X  +  c(t)X2  (19) 


Eqs.  (1),  (2),  (4)  and  (7)  can  be  used  to  determine  the  fol¬ 
lowing  equation  for  the  overpotential  in  dimensionless  form: 


2 


dr)  dr) 
~dr  ~  ax2 


2  ^ 

v  rj 


(id 


where  v2  (dimensionless  exchange  current  density)  [6]  is 


_  aioiaa  +  ac)FLz  ^11 


RT 


a 


K 


(12) 


and 


T) 


* 


r)F 

RT 


t 

aC(\\  (1  /k  +  l/a)L2’ 


(13) 


The  corresponding  dimensionless  initial  and  boundary  con¬ 
ditions  are 


rf  =  0  at  r  =  0  and  for  0  <  X  <  1 

(14) 

=  8  at  X  =  0  and  for  r  >  0 
dA 

(15) 

=  —8/3  at  X  =  1  and  for  r  >  0 

dx  F 

(16) 

where  f3  is  the  ratio  of  the  solution  phase  to  matrix  phase 
conductivity  (/3  =  k/o)  and  8  is  the  dimensionless  current 

density: 

s  =  i(FL) 

V  kRT  ) 

(17) 

2.7.  Exact  transient  solution 


where  a(r),  b(x)  and  c(r)  may  be  functions  of  the  dimen¬ 
sionless  time.  Two  of  these  time-dependent  parameters  (b(r) 
and  c(r))  can  be  determined  using  the  boundary  conditions. 
By  applying  Eq.  (19)  to  the  boundary  condition  at  X  =  0 
we  get 

b(  r)  =  5  (20) 


Similarly,  c(r)  can  be  obtained  by  using  the  boundary  con¬ 
dition  at  X  =  1  as  follows: 

c(  T)  =  -^(l+/J)  (21) 


An  average  77*  can  be  determined  by  using  Eq.  (19): 


r) 


* 

avg 


dX  =  a(x)  + 


Hr) 

2 


+ 


c(r) 

3 


(22) 


Eq.  (11)  together  with  Eqs.  (19)  and  (21)  can  be  used  to 
derive  an  equation  for  the  time  dependence  of  the  dimen¬ 
sionless  average  overpotential,  77*vg.  First,  use  Eqs.  (19)  and 
(21)  to  obtain  an  expression  for  the  second  derivative  of  77* 
with  respect  to  X.  Next,  put  this  expression  into  Eq.  (11)  and 
integrate  both  sides  of  the  resulting  equation  with  respect  to 
X  from  X  =  0  to  1 .  The  result  is 


d  77* 
'a'/ avg 

dr 


-8{\+P)-v2i1l/% 


(23) 


Eq.  (23)  can  be  solved  by  using  the  initial  condition  77*  =  0 
at  r  =  0  for  all  X  to  obtain 


*  _  _v2T8(l+/3)  8( l+P) 

7?avg  e 


V" 


(24) 


Next,  the  time-dependent  parameter  a( r)  can  be  determined 
by  substituting  Eqs.  (20),  (21)  and  (24)  into  Eq.  (22).  The 
result  is 


,  ,  _v2t8{\+P)  8{l  +  P)&p  8 

a(x)  =  e  T - ■= - ■= - h  —  -  - 


(25) 


An  exact  transient  solution  for  this  model  can  be  derived 
using  the  separation  of  variables  method  [7]  and  is  given  by 

*  _  5(1  +  fi)  e“l’2'r  5[cosh(v[l  -  X])  +  ^cosh(vX)] 

^  v2  v  sinh(v) 

00 

+2 Sj^An  cos (nitX)  c~(nV+v2)T  (18) 

n  —  1 

where  An  =  (/3cos(nn)  +  l)/(v2  +  n2 ix1) 

The  number  of  terms  required  for  convergence  of  this  se¬ 
ries  can  be  significant  and  depends  strongly  on  the  dimen¬ 
sionless  current  density  8. 


The  dimensionless  overpotential  based  on  the  polynomial 
approximation  (simplified  dynamic  model)  is  obtained  by 
substituting  a( r),  b(x)  and  c(x)  into  Eq.  (19).  The  result  is 


77  = 


^  8/3  _  2  8(1  +  p) 

- h  —  +e  y  T - 


5(1  +  P) 


vx 


vx 


+Sx-8X±FX2 


(26) 


Comparing  Eq.  (18)  to  Eq.  (27)  shows  that  the  simpli¬ 
fied  dynamic  model  does  not  contain  an  infinite  series  and 
is,  consequently,  much  simpler.  Note  that  the  simplified  dy¬ 
namic  model  involves  all  the  physical  and  kinetic  parameters 
of  the  model. 
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2.3.  Potential  drop  across  the  electrode 


transient  solution  (Eq.  (18))  as 


Combining  Eqs.  (1),  (4)  and  (7)  gives  the  equation  for 
solid  phase  potential  as 

d2<pi 


dr]  aio(aa+ac)F 

'7a?=‘,c‘ll¥  + — rt — " 


(27) 


Substitution  of  the  dimensionless  variables  as  defined  in 
Eq.  (13)  into  Eq.  (27)  and  letting 

<PlF 


<t>\  = 


RT 


(28) 


yields  an  equation  for  the  dimensionless  solid  phase  potential 

)2  a* 

n  an  r,  . 

(29) 


d2<t> * 

dX2 


1  + 


3  ri*  2  * 

+  V1J 


dr 


Using  Eq.  (11),  Eq.  (29)  can  be  rewritten  as 

12  ^  p  d2i]* 


fn  = _ 

dX2  l+fi  dX2 


(30) 


The  boundary  conditions  for  <fi*x  and  77*  in  Eq.  (30)  are 


* 


#1 

dx 

90* 


=  0  and 


di] 


* 


dx 


=  S  at  X  =  0 


=  —8/3  and 


dr] 


* 


=  -Sp  at  X  =  1 


(31) 


(32) 


dx  dx 

Eq.  (30)  can  be  integrated  once  with  respect  to  X  to  yield 
90*  p  dr]* 


dX 


1  +  p  dX 


+  C 


(33) 


where  C  is  an  integration  constant  that  can  be  determined 
using  one  of  the  boundary  conditions  (Eq.  (31)  or  (32)): 

sp 


c  =  - 


1  + 


(34) 


Integrating  Eq.  (34)  between  the  limits  X  =  0  to  1  gives  the 
change  in  the  dimensionless  solid  phase  potential  across  the 
electrode 


0*  I X=1  —  01  lx=0  = 


(r]*\x=i-r]*\x= 0)- 


* 


1  + 


Sp 


1  + 


(35) 


The  potential  drop  across  the  porous  electrode  can  be 
defined  as 


V*  =  (tf|x=i-$lx=o) 


(36) 


.*  8(1+ p)c 


—  V2T 


y*  = 


sp 


S[(l  +  p2)  cosh(v)  +  2/3] 


V" 

00 


+29^  Bn 


1  + 


,  — (ft2  7T2+V2)t 


v  sinh(v)(l  +  P) 


(39) 


a = 1 


where 


Bn  = 


( Pcos(nTt )  +  1)^ 
v 2  +  n2it2 


The  dimensionless  voltage  drop  across  the  porous  elec¬ 
trode  from  X  =  1  to  X  =  0  as  determined  by  using  the 
simplified  dynamic  model  (Eq.  (27))  consists  of  only  two 
terms: 


,*  9(l+£)[e-^-l]  5(1 +j8) 


y^  = 


vx 


(40) 


2.4.  Interfacial  current  density 


Substituting  for  i\  in  Eq.  (4),  from  Eq.  (1)  gives  the  in¬ 
terfacial  current  density  as 


]n  — 


1  d2<t>i 

dx2 


aa 


(41) 


Introducing  the  dimensionless  variable  0]*  (Eq.  (28)  and  X 
Eq.  (13))  into  Eq.  (41)  gives  the  dimensionless  interfacial 
current  density  (j*).  j*  can  be  written  in  terms  of  the  di¬ 
mensionless  overpotential  using  Eq.  (30): 


2 


/*  = 
Jn 


1  dAr] 


P  +  1  dX2 


(42) 


Evaluating  the  second  derivative  of  the  dimensionless  over¬ 
potential  with  respect  to  the  distance  from  Eq.  (18)  gives  the 
expression  for  j*  as 


•* 


Jn 


EP  COS(t27T) H- 1  9  9  /„2_2_i_1,2\_ 

—4 — 2  -  ;  - n2jt2  cos(uTrX)  e“(n  71  +v  )r 

(v2+n2n2)(P+\) 


ft=  1 


Sv  [cosh(v[l  —  X])  +  p  cosh(v2Q] 
(P  +  1)  sinh(v) 


(43) 


3.  Results  and  discussion 


or 


y*  =  (0*|x=!  -  01  lx=o)  +  9*lx=0 


(37) 


By  substituting  Eq.  (35)  into  Eq.  (37),  we  get  the  dimen¬ 
sionless  voltage  drop  across  the  electrode  as 


0*\x=o  +  Po*  \x=i  ~  jg 

1  +  P 


(38) 


This  dimensionless  voltage  across  the  porous  electrode  (y*) 
from  X  =  1  to  X  =  0  can  be  determined  by  using  the  exact 


From  the  equations  presented  above,  it  is  clear  that  both 
overpotential  and  voltage  drop  across  the  porous  electrode 
are  functions  of  the  dimensionless  current  density  9,  di¬ 
mensionless  exchange  current  density  v2  and  the  ratio  of 
conductivities  p.  The  distribution  of  the  overpotential  across 
the  porous  electrode  obtained  from  the  transient  solution 
is  compared  to  the  solution  obtained  using  the  simplified 
dynamic  model  in  Fig.  2  for  the  parameter  values  9  =  —  1, 
v2  =  P  =  1.  This  figure  also  shows  the  overpotential  dis¬ 
tribution  at  different  dimensionless  times.  The  change  in 
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Fig.  2.  Distribution  of  over  potential  across  the  porous  electrode 
(v2  =  ft  =  1,  8  =  —  1).  The  exact  solution  (continuous  line)  and  the 
simplified  dynamic  model  solution  (symbols)  agree  well. 


voltage  drop  across  the  porous  electrode  with  time  is  plotted 
using  both  the  exact  solution  and  the  simplified  dynamic 
solution  in  Fig.  3.  The  voltage  increases  linearly  at  short 
times  and  levels  off.  The  voltage  profiles  are  similar  to 
the  transient  potential  response  of  insertion  electrodes  with 
linear  kinetics  and  double-layer  charging  as  shown  by  Ong 
and  Newman  (Fig.  5  of  ref.  [5]). 

It  is  worth  noting  that  Eq.  (18)  can  be  simplified 
as  Eq.  (44),  when  the  pseudocapacitor  involves  only 


T 

Fig.  3.  Change  in  voltage  across  the  porous  electrode  with  dimensionless 
time  (v2  =  /3  =  1).  The  exact  solution  (continuous  line)  and  the  simplified 
dynamic  model  solution  (symbols)  agree  well. 


double-layer  charging  (i.e.,  vx 


T]  = 


8  o  8/3 

-5(1  +  P)t  -  ~(3X2  -6X  +  2)-  —{3X‘ 

6  6 

°°  cos(rutX)(/3  cos  (rut)  +  \)c~n2jr2T 


-1) 


n= 1 


n2n2 


(44) 


It  can  be  shown  that  Eq.  (44)  is  equivalent  to  Eq.  (16) 
given  by  Srinivasan  and  Weidner  [2]  provided  the  bound¬ 
ary  conditions  and  the  dimensionless  variables  given  here 
are  redefined  according  to  those  given  in  ref.  [2].  Similarly 
when  v2  — >►  0,  Eq.  (39)  reduces  to  the  voltage  expression 
given  in  Eq.  (45): 

y*  =-8r  G8+D  — ^(j8+l) 

~  250Scos(«7r)  +  l)2  Vr 

nWtf+l) 

n= 1 


Eq.  (45)  is  equivalent  to  Eq.  (15)  of  ref.  [2]  when  the  initial 
voltage  on  discharge  is  assumed  to  be  equal  to  zero.  Hence 
the  transient  model  developed  for  the  pseudocapacitor  re¬ 
duces  to  that  of  the  double-layer  capacitor  when  the  kinetics 
of  the  faradaic  reaction  is  negligible  (v2  — >  0).  The  voltage 
profiles  for  a  pseudocapacitor  are  compared  to  that  of  a 
double-layer  capacitor  for  different  values  of  the  discharge 
current  density  in  Fig.  4.  Newman  pointed  out  that  for  large 
values  of  the  dimensionless  exchange  current  density  (v2) 
and  the  dimensionless  applied  current  density  (5),  the  ohmic 
effects  dominate  and  the  reaction  distribution  is  non-uniform 


T 

Fig.  4.  The  effect  of  dimensionless  exchange  current  (v2)  density  on  the 
voltage  profile  across  the  porous  electrode  for  /3  =  0.  The  dimensionless 
applied  current  density  (5)  are  as  indicated  in  the  figure.  Dotted  lines 
are  the  voltage  profile  for  v2  =  0.01  simulated  using  the  exact  solution. 
Solid  lines  are  the  voltage  profile  for  v2  =  1  simulated  using  the  exact 
solution.  The  exact  solution  (line)  and  the  simplified  dynamic  model 
solution  (symbols)  agree  well. 
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Fig.  5.  Distribution  of  dimensionless  interfacial  current  density  (j*)  for 
v2  =  /3  =  l,  8  =  —  1  and  at  different  dimensionless  times,  j*  (Eq.  (43)) 
is  derived  using  the  exact  solution  for  the  overpotential. 


Fig.  7.  The  effect  of  dimensionless  applied  current  density  on  the  distri¬ 
bution  of  dimensionless  interfacial  current  density  (j*)  for  v2  =  /3  =  1, 
8  =  —  1.  j*  (Eq.  (43))  is  derived  using  the  exact  solution  for  the  overpo¬ 
tential. 


[6].  The  dimensionless  interfacial  current  (j*)  distribution 
for  v2  =  1  is  illustrated  in  Fig.  5  at  different  r.  The  distribu¬ 
tion  of  j*  reaches  a  steady  state  when  r  >  0.15.  Compared 
to  the  steady  state  j*  distribution  for  v2  =  0.01  (small  value 
of  the  dimensionless  exchange  current  density)  as  shown  in 
Fig.  6,  it  is  clear  that  the  steady  state  reaction  distribution  is 
non-uniform  for  v2  =  1 .  The  steady  state  distribution  of  j* 
for  v2  =  0.01  (Fig.  6)  is  more  uniform  and  corresponds  to 
reaction  distribution  of  the  double-layer  capacitors  [2].  The 
effect  of  8  on  the  reaction  distribution  is  given  in  Fig.  7.  For 


higher  values  of  the  discharge  current  density  the  reaction 
distribution  becomes  non-uniform.  The  effect  of  the  conduc¬ 
tivity  ratio  (P)  on  the  reaction  current  density  is  to  shift  the 
reaction  from  one  face  to  another  [6].  This  is  evident  from  the 
plot  of  the  reaction  distribution  for  low  value  of  ft  as  given 
in  Fig.  8.  For  ft  =  0,  solid  phase  conductivity  is  higher  than 
the  solution  phase  conductivity,  the  reaction  rate  is  higher 
at  the  solution/porous  electrode  interface.  The  reaction  dis¬ 
tribution  for  p  =  10,  where  the  solution  phase  conductivity 
is  higher  than  the  solid  phase  conductivity  the  profiles  are 


Fig.  6.  Distribution  of  dimensionless  interfacial  current  density  (j*)  for 
v 2  =  0.01,  f  =  1,  8  =  —  1  and  at  different  dimensionless  times,  j* 
(Eq.  (43))  is  derived  using  the  exact  solution  for  the  overpotential. 


Fig.  8.  Distribution  of  dimensionless  interfacial  current  density  (j*)  for 
v2  =  1,  /3  =  0,  8  =  —  1  and  at  different  dimensionless  times,  j*  (Eq.  (43)) 
is  derived  using  the  exact  solution  for  the  overpotential. 
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a  mirror  image  of  the  plots  for  /3  =  0.  Hence  the  reac¬ 
tion  rate  is  higher  at  the  porous  electrode/current  collector 
interface. 

4.  Conclusion 

An  analytical  solution  for  a  macroscopic  model  of  a 
porous  electrode  with  linear  kinetics  was  compared  to  a 
simplified  dynamic  model.  Comparing  the  results,  it  was 
shown  that  the  simplified  dynamic  or  approximate  model 
agrees  well  with  the  exact  transient  model  for  the  cho¬ 
sen  parameters  (v2  =  /3  =  1).  The  model  includes  both 
double-layer  charging  and  a  faradaic  process  approximated 
by  linear  kinetics.  The  exact  transient  analytical  solution  for 
the  overpotential  was  used  to  determine  the  dimensionless 
interfacial  current  density  (Eq.  (43)).  The  dimensionless 
interfacial  current  density  distribution  is  analyzed  as  a 
function  of  the  parameters  (v2,  /3  and  8).  A  rigorous  error 
analysis  based  on  the  parameters  8 ,  /3  and  v  could  be  done 
as  in  ref.  [7].  This  analysis  would  determine  the  range  of 
the  parameters  for  which  Eq.  (27)  could  be  used  safely.  In 


addition,  polynomial  approximations  of  higher  order  [8] 
could  be  used  to  increase  the  accuracy. 

Acknowledgements 

The  authors  are  grateful  for  the  financial  support  of  the 
project  by  National  Reconnaissance  Organization  (NRO)  un¬ 
der  contract  #  NRO  000-03-C-0122. 


References 

[1]  R.  Kotz,  M.  Carlen,  J.  Electrochem.  Soc.  45  (2000)  2483. 

[2]  V.  Srinivasan,  J.W.  Weidner,  J.  Electrochem.  Soc.  146  (1999)  1650. 

[3]  C.  Lin,  J.A.  Ritter,  B.N.  Popov,  R.E.  White,  J.  Electrochem.  Soc.  146 
(1999)  3168. 

[4]  D.  Dunn,  J.  Newman,  J.  Electrochem.  Soc.  147  (2000)  820. 

[5]  I.J.  Ong,  J.  Newman,  J.  Electrochem.  Soc.  146  (1999)  4360. 

[6]  J.S.  Newman,  Electrochemical  Systems,  2nd  ed.,  Prentice-Hall,  En¬ 
glewood  Cliffs,  NJ,  1991. 

[7]  V.R.  Subramanian,  R.E.  White,  J.  Power  Sources  96  (2001)  385. 

[8]  V.R.  Subramanian,  J.A.  Ritter,  R.E.  White,  J.  Electrochem.  Soc.  148 
(2001)  E444. 


